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We present first-principles calculations to predict several three dimensional (3D) topological in- 
sulators in quaternary chalcogenide compounds which are made of I2-II-IV-VI4 compositions and 
in ternary compositions of I3-V-VI4 famatinite compounds. Among the large members of these two 
families, we give examples of naturally occurring compounds which are mainly Cu-based chalco- 
genides. We show that these materials are candidates of 3D topological insulators or can be tuned 
to obtain topological phase transition by manipulating the atomic number of the other cation and 
anion elements. A band inversion can occur at a single point V with considerably large inversion 
strength, in addition to the opening of a bulk band gap throughout the Brillouin zone. We also 
demonstrate that both of these families are related to each other by cross-substitutions of cations 
in the underlying tetragonal structure and that one can suitably tune their topological properties 
in a desired manner. 
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I. INTRODUCTION 

The recent discovery of topological insulators 1 - - — has realized a long-sought opportunity of implementing exotic 
quantum phenomenon in practical materials J^r— However, existing topological materials are restricted to binary 
and ternary compositions of heavy elements^— *L2r— , which limits the possibility of combining correlated electronic, 
magnetic, superconducting and other local order properties with topological phenomenon. It is highly desirable, 
therefore, to expand the menu of available materials into the domain of compounds which are relevant for device 
fabrication with tuning capabilities and accessible crystalline cleavage. 

Topological insulator order and local order of spontaneous symmetry breaking such as magnetism and supercon- 
ductivity has recently been classified in different spatial dimensions and with different discrete symmetries. The band 
topology in the two-dimensional quantum spin-Hall effect or three dimensional non-trivial topological insulators are 
required to possess a time-reversal invariant number Z2 = -I. Materials with strong spin-orbit coupling can host 
non-trivial topological phases in both two and three spatial dimensions when the time-reversal symmetry remains 
invariant. The topological insulating materials realized to date are based on the binary and ternary combinations 
of mostly heavy metals ranging from HgTe^ Bi/Sb-based binary systems£ii2~— to Tl-based ternary compounds 28 . 
However, the principle allure of topological insulators stems from their potential technological applications, which 
require the time-reversal-invariant topological ground state to be placed along with the broken symmetry of local 
orders* 1 ^ - — Of particular interest are the magnetic atom doped topological insulators for realizing effects of magnetic 
impurities and ferromagnetism on the topological surface states.— ~— Furthermore, the realization of elementary ex- 
citations that satisfy non-Abelian statistics — the so-called Majorana fermions, important for topological quantum 
computers, requires proximity of superconducting and ferromagnetic insulating phasesj 15 ' 32 i 33 Therefore, it is impor- 
tant to search for material classes which provide more chemical and structural flexibility to tune the ground state 
properties using cross-substitution of elements. 

Here, we report the theoretical prediction of three dimensional topological insulator in the pristine state of ex- 
isting materials in quaternary chalcogenide compounds which are made of I2-II-IV-VI4 compositions and in ternary 
compositions of I3-V-VI4 famatinite compounds. Both families are related by cross-substitutions of cations in the 
underlying tetragonal structure. Quaternary semiconductors, especially, exhibit more flexible properties arising from 
their enhanced chemical and structural freedom and provide a handle on the surface quantum control through the 
large possibility of substitution with the magnetic, non-magnetic and other correlated electronic elements. Through 
the ternary-to-quaternary mutation and control of the atomic configuration, one can perform band engineering to 
tailor the topological quantum phenomenon for transport and quantum information processing application. 

We begin the discussion with the ternary composition of I3-V-VI4 compounds which are known to crystallize in 
famanitine structures and are basically of the sphalerite-type. These materials belong to the space group /42m, in 
which the group VI atoms are surrounded by three group I and one group V elements. As a result, they obey the 
octet rule and form a (I-VI)3(V-VI) superlattice structure, see Fig. 1(a) for CusSbSe4 as an example. In contrast 
to the high-symmetry zinc-blende sublattice of ternary compound s 28 ' 34 ! 35 , famatinite structure naturally achieves a 
tetragonal lattice distortion along the c— axis (c < 2a) due to the turning on of the strong interlayer coupling between 
the two cation planes. This also results in a mismatch between the cation-anion bond lengths in two zinc-blende 
formula units of the unit cell which helps lower the total energy of the famatinite ground state phase. 

Famatinite compounds I3-V-VI4 evolve to the quaternary chalcogenides I2-II-IV-VI4, when one of the group-I 
element of the former compounds is replaced by one group-II element and the group-V element is changed to group- 
IV element; see Fig. 1(b) for example Cu2CdSbSe4. High-resolution transmission electron microscopy and X-ray 
diffraction analysis shows that these materials have tetragonal crystallographic structure which belongs to /42m 
space group. The structure can be described by a (I-II) 2 (II-VI)(IV-VI) sublattice of the two zinc-blende formula units 
as before^ The substitution of larger atoms in the tetragonal basis of these systems helps expand the lattice and 
thus both the tetragonal distortion and the crystal-field splitting increases considerably. Although these changes are 
relatively small, nevertheless, quaternary chalcogenides have more structural freedom and more complicated electronic 
and chemical properties than their binary or ternary counterparts. 

The typical topological phases can be classified as trivial band insulator, a non-trivial topological semimetal (zero- 
bandgap) and a non-trivial topological insulator(finite-bandgap). The corresponding band symmetries near T point 
are expressed schematically in Fig. lc-e respectively. For a trivial topological phase with the Z2 topological order 
equal to +1 in a zinc-blende crystal such as CdTe, a bandgap exists at the Fermi-level as well as the twofold degenerate 
s— like states marked by darked-blue dot lies in the conduction band. In a non-trivial topological phase case such 
as HgTe, band inversion occurs whereas s— like states drop below the fourfold degenerate p— like states of J = 3/2 
multiplet, with Z2—-I. The fourfold degeneracy is due to the cubic symmetry in the zinc-blende structure. To 
achieve a topological insulating state, the degeneracy is required to be lifted and a band gap is opened (Fig. 1(e)). 
Since the non-cubic famatinites and quaternary chalcogenides are derived from zinc-blende structures, the non-trivial 
topological insulating phase could be realized in these materials. Indeed, our first-principles calculations predict 
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FIG. 1. Crystal structure and the property of the topological band-inversion, (a), The crystal structure of ternary 
famatinite compounds is given for a representative material Cu3SbSe4. (b) The crystal structure of quaternary chalcogenide 
compounds is illustrated for Cu2CdSbSe4 as an example, (c)-(e) Band structures illustrated schematically near the time-reversal 
invariant point V for a trivial band insulator in (c), a non-trivial topological semimetal in (d) and a non-trivial topological 
insulator in (e). Dark-blue dots represent the twofold degenerate s— like states. In the non-trivial case, band inversion occurs 
when s-like states drops below the four-fold degenerate p— like states (J = 3/2) (blue lines). Owing to the lattice distortion, 
the degeneracy of the J = 3/2 states vanishes in the present non-cubic crystal structures. 

that 3D non-trivial topological insulating phase exists in the compounds of famatinite and quaternary chalcogenides 
families. 



II. METHODS 

The crystal structures of ternary famatinite compounds and quaternary chalcogenide compounds are taken from 
literature^ - — . First-principles band calculations were performed with the linear augmented-plane-wave (LAPW) 
method using the WIEN2K package^ in the framework of density functional theory (DFT). The generalized gradient 
approximation (GGA) of Perdew, Burke, and Ernzerhof 5 - was used to describe the exchange-correlation potential. 
Spin orbital coupling (SOC) was included as a second variational step using a basis of scalar-relativistic eigenfunctions. 



III. RESULTS 

The bulk band structures along high symmetry path M(tt,it) — r(0, 0) — X(ir,0) are shown in Fig. 2 for two 
representative compunds Cu3SbS4 and Cu3SbSe4 which belong to the famatinite family. The effect of the zone- 
folding is clearly evident at the high-symmetry point T. In the ground state of these compounds, the structural 
compression along c-axis (i.e. c < 2a) compared to cubic zinc-blende lattice lifts the four-fold degeneracy of the 
cation p— states at the T— point owing to the crystal-field splitting. In Cu3SbS4, an insulating energy gap between 
the p— states is present at the Fermi level throughout the Brillouin zone as shown in Fig. 2(a). However, as the 
nuclear charge of the anion is increased from S atom to Se atom, such substitution causes the conduction band to 
drop below the Fermi level at the X-point, yielding an electron pocket, see Fig. 2(b). Simultaneously, the valence 
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FIG. 2. The bulk band dispersion of famainite compounds with non-trivial topological phases, (a)-(b) Electronic 
structure of Cu3SbS4 and Cu3SbSe4. In both figures, the black dots are the p— states of the cation sites while the size of the blue 
dots is proportional to the probability of s— orbital occupation on the anion site. The gray shadings highlight the low-energy 
region near T— point where the inverted band-inversion has occurred. 



band maximum also gradually moves above the Fermi-level making Cu3SbSe4 a metal. Nevertheless, a finite direct 
gap persists throughout the Brillouin zone and thus the Zi topological invariant can still be defined for the valence 
bands in these two materials and the inverted band order is retained at only T— points. Therefore, the Cu3SbS4 and 
Cu3SbSe4 are topologically non-trivial insulator and metal, respectively. 

The mutation of the crystal structure from the ternary to a quaternary compound is also evident in the electronic 
structure of these two families of materials as the overall changes in the electronic structure are only subtle near the 
energy region of present interests. Due to the tetragonal symmetry, all the Cu-based quaternary materials listed in 
Fig. 3 host a gap in the p— states close to the Fermi-level, where the gap magnitude depends on the extra cation 
cross-substitutions. With increasing atomic mass from Zn to Cd to Hg atom, this gap decreases mainly due to the 
downshift of the s— like conduction bands at T— point. The Fermi-level goes through the energy gap of the p— states 
in most materials, harboring a topological insulating state in these compounds. However, only in two compounds 
Cu2ZnGeTe4 and Cu2CdGeTe4, the low-energy electronic properties reveal that the conduction band minimum is 
shifted to the X point and drops below the Fermi level to form bulk electron pockets. Simultaneously, the concave- 
upward shaped valence band maximum pops up above the Fermi level at T— point, making these systems intrinsically 
bulk metallic. We find very similar nature of the band-inversion across these materials in which the split p— states lie 
in energy above the twofold-degenerate s— states at a single time-reversal momentum T— point, representing a band- 
inversion relative to the natural order of s— and p— type orbital derived band structure. Therefore, these materials are 
intrinsically bulk non-trivial topological insulators or semimetals, except Cu2ZnGcS4, Cu2ZnSnS4, and Cu2HgGeS4 
which are trivial topological insulators. 

Finally, we demonstrate our tuning procedure to generate a topological phase transition from a trivial band insulator 
to the non-trivial topological insulator. In the present study, we take advantage of the structural freedom of the solids 
in ternary famatinite and quaternary chalcogenides to achieve non-trivial topological phases. In the Fig. 4, we discuss 
one representative example of this route here in going from the quaternary chalcogenide Cu2ZnSnS4 as the starting 
point to the ternary famatinite CusSbS4 compound. Note that the topological phase transition is a generic feature and 
there exists numerous such route. Cu2ZnSnS4 is a trivial band insulator whereas Cu3SbS4 is a topological insulator 
in their pristine conditions as discussed before in Figs. 3(i) and 2(b), respectively. We systematically manipulate the 
overall crystal structure of our starting compound by making the atomic number Z of the constituent atoms, the 
lattice constant and the internal displacement of the anion as variable parameters. 

To illustrate this process, we denote the atomic number Z of various elements by Zmi, Zm2, Zmz, and Zma- For 
our starting element Cu2ZnSnS4, we have Zmi=29, Zm2=30-%> Zm3=50+x, Zm4=16, where x varies between to 1 
in the phase transformation process. At each step of changing the parameter x, we also pay attention to the lattice 
constant and the internal displacements of the anion by keeping the same ratio of changes corresponding to the change 
in x of atomic number Zs from Cu2ZnSnS4 to Cu3SbS4. Note that variation of the atomic number Zs can also be 
considered as a doping effect in the virtual crystal approximation while the total nuclear charge is neutral during the 
whole process. 

We monitor the band structure along high symmetry path M(ir, 7t) —> T(0, 0) — > X(ir, 0) at every step of the phase 
transition process from Figs. 4(a)-(e). While the overall band- structure remains very much the same throughout the 
process, the band characters and the band gap vary dramatically. In Fig. 4(b), on increasing x by 0.25, the size of 
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FIG. 3. Electronic structure of quaternary chalcogenide compounds. Bulk electronic structure of several Cu-based 
quaternary chalcogenides compounds as indicated in the title of each figures. The black dots, blue dots and the gray shading 
have the same meaning as in Fig. 2. The compounds without a gray-shaded area are trivial topological insulators without any 
inverted band symmetries whereas the others are non-trivial topological insulators or semimetals with Z-z = — 1 topological 
order. 
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FIG. 4. Topological phase transition between quaternary chalcogenide and famainite compounds. Bulk electronic 
structure of the compounds at various stages of the topological phase transition process by tuning the atomic number Z, the 
lattice constant and the internal displacement of the anion are shown here, (a) and (b) are the trivial band insulator, (c) the 
critical point of the topological phase transition whereas (d) and (e) are the non-trivial topological insulators. 



band gap at T starts to shrink compared to the parent compound in Fig. 4(a), although the s— character continues to 
dominate in the conduction band (blue dots) . The critical point of the phase transition is reached when the value of x 
is increased to 0.5 in Fig. 4(c), at which point the bottom of the conduction band and the top of the valence band touch 
at a single momentum point at T leading to a critical topological point. A full band inversion is achieved in Fig. 4(d) 
which recreates an insulating band gap but topologically non-trivial in nature and the s— like orbital character shifts 
much below the valence band. Finally, in Fig. 4(e), when x equals 1, Cu3SbS4 is a non-trivial topological insulator. 



IV. CONCLUSION 

We have shown via first-principle calculations that the ternary famatinite and quaternary chalogenide compounds 
are non-trivial topological insulators in their pristine phase. The involvement of partially filled and /— electrons 
in these compounds naturally opens up the possibility of incorporating magnetism and superconductivity within 
the topological order. The large-scale tunability available in the quaternary compounds adds versatility in using 
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topolgical insulators in multifunctional spin-polarized quantum and optical information procesing applications. No- 
tably, the copper-based quaternary chalcogenides are known to have non-linear optoelectronics and thermodynamics 
applications.—. 

During the final stages of our work, we became aware of the work of Su-Huai Wei et oL— who have also proposed 
that some of the multinary chalcogenides are potential candidates to be non-trivial topological insulators. 
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